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; ABSTRACT 

, Deep 1.1 mm continuum observations of 1E0657-56 (the "Bullet Cluster") taken with 

C*~) ' the millimeter-wavelength camera AzTEC on the 10-m Atacama Submillimeter Tele- 

| scope Experiment (ASTE), have revealed an extremely bright (Si.imm = 15.9 mJy) 

unresolved source. This source, MMJ065837-5557.0, lies close to a maximum in the 
density of underlying mass-distribution, toward s the larger o f the t wo interacting clus- 
ters as traced by the weak-lensing analysis of IClowe et al.l ((2006) . Using optical-IR 
colours we argue that MMJ065837-5557.0 lies at a redshift of z = 2.7 ± 0.2. A lensing- 
derived mass-model for the Bullet Cluster shows a critical-line (caustic) of magnifi- 
cation within a few arcsecs of the AzTEC source, sufficient to amplify the intrinsic 
millimeter-wavelength flux of the AzTEC galaxy by a factor of ^> 20. After subtrac- 
tion of the foreground cluster emission at 1.1mm due to the Sunyaev-Zel'dovich effect, 
and correcting for the magnification, the rest-frame FIR luminosity of MMJ065837- 
5557.0 is ^5 10 12 L Q , characteristic of a luminous infrared galaxy (LIRG). We explore 
various scenarios to explain the colors, morphologies and positional offsets between 
the potential optical and IR counterparts, and their relationship with MMJ065837- 
5557.0. Until higher-resolution and more sensitive (sub)millimctcr observations are 
available, the detection of background galaxies close to the caustics of massive lensing 
clusters offers the only opportunity to study this intrinsically faint millimeter-galaxy 
population. 
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1 INTRODUCTION 

A high-redshift (z ^> 1) strongly-evolving submillimeter 
galaxy population (hereafter SMGs) has been identified 
in a series of continuum imaging surveys during the last 
decade. The most significant s urveys have been u nder- 
taken with SCUBA at 850^im (jCoppin et all (|2006h and 



references therein), M AMBO at 1.2m m (iGreve et all 12004 ; 
iBertoldi et all 12007ft . BOLOCAM (jLaurent et all 12005 ) 
at 1.1 mm, and more recently with AzTEC at 1.1mm 
(|Scott et alj|200§ ). These collective submillimeter surveys, 
together with comprehensive multi-wavelength imaging and 
spectroscopic follow-up spanning X-ray to radio wave- 
lengths, have demonstrated that SMGs are representative 
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of a population of massive, dust-enshrouded, optically-faint 
galaxies undergoing significant starformation with rates ^> 
200M Q /yr _1 . The median redshift of the bright SMG popu- 
lation, based on spectroscopic and photometric redshifts, is 
~ 2.5, w ith approximately 50% of the population at 1.9 < 
z < 2.9 (IChapman et alj|2003t |2005| ; lAretxaga et~ail 12001 



120071 : iPope et al.ll2005h . Scaling the observed flux-densities 
of all SMGs with detections in the above flux-limited sur- 
veys to 850/xm, the majority of the population have been 
detected with flux densities of 2 mjy < Ssso fan < 12 mjy. 

Making the connection between ultra-luminous SMGs 
and their local analogs - presumably UL1RGS and L1RGS 
- requires either larger telescopes such as the Large Mil- 
limeter Telescope (LMT) and the Atacama Large Millime- 
ter Array (ALMA), or using the presence of foreground 
structure to amplify the faint background galaxies. To date, 
the deepest SMG surveys have been conducted towards 
lensing-clusters that amplify an intr insically fainter popula- 



tion of SMGs in the high-z universe (ISmail et al.lll997l 
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the h i gh- z u niverse (binail et al. II 199 /I. 
20021; ICowie et al.ll2002l; iKneib et al.1 
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2004 



Knuds en et al . 2006). Imaging massive clusters offers the ad- 



vantages of amplification of the background universe due 
to the gravitational-potential of the cluster and, in general, 
an increased effective resolution in the source plane (and 
correspondingly decreased survey confusion limit). Unfortu- 
nately, this same advantage for those studying faint back- 
ground galaxies comes at a cost to those interested in study- 
ing the clusters themselves using the Sunyaev-Zel'dovich Ef- 
fect (SZE). The amplified (sub)millimeter- wavelength emis- 
sion from members of the SMG field population that are 
lensed by the large-scale cluster potential and by the in- 
creased probability of galaxy-galaxy lensing along the line 
of sight to the cluster, increase the likelihood of point-source 
contamination of the SZE. 

In this paper we present the discovery of an extremely 
bright millimeter-selected SMG (MMJ065837-5557.0) in 
the direction of the massive X-ray luminous Bullet Clus- 
ter (z=0.297). The 1.1 mm-wavele ngth observations w ere 
made with the AzTEC instrument l| Wilson et al.l 120081 ) on 
the Atacama Subm illimeter Telescope Experiment (ASTE, 
lEzawa et all l|2004h ). These AzTEC data were undertaken 
as part of a larger project to trace the evolution in the 
surface-density of SMGs towards over-dense regions in the 
low and high-redshift Universe. Using archival HST and 
Spitzer data, we identify potential optical and IR counter- 
parts to the AzTEC source. The photometric-colors, and the 
small positional-offset of the AzTEC source and the optical- 
IR counterparts from the critical magnification-line derived 
from cluster mass-models strongly support the suggestion 
that MMJ065837-5557.0 is a moderately-luminous high-z 
galaxy that has been lensed and amplified by the massive 
foreground Bullet Cluster. 

Throughout this paper we adopt the following cosmo- 
logical model: a Hubble constant Ho = 72 kms - Mpc~ , 
and density parameters fl m = 0.3 and J1a = 0.7. 



2 AZTEC/ ASTE OBSERVATIONS 

AzTEC is a 144-element bolometric receiver cu rrently tuned 
to op erate in the 1.1mm atmospheric window (|Wilson et al.l 
2008). We have previously used AzTEC to complete a suc- 



cessful set of observations of the submillimeter galaxy pop- 
ulation in blank-fields fro m the 15 m d iamet er James Clerk 
Maxw ell Telescope (e.g. IScott et all (|2008l ); iPerera et ail 
(2008); Austermann et al, in preparation). In 2007 AzTEC 
was mounted on the ASTE telescope, a 10 m diameter an- 
tenna located on Pampa La Bola, near Cerro Chajnantor 
in Chile. ASTE provides AzTEC with an angular resolution 
of 30" FWHM. One season of observations (June-October 
2007) from this excellent high (~ 4900m) and dry site in the 
Atacama Desert has been completed and another season will 
begin in July 2008. 

The bright millimeter- wavelength source, shown in Fig- 
ure [T] was found in our 202 sq. arcmin AzTEC imaging sur- 
vey of the Bullet Cluster and its environment (Ezawa et al., 
in preparation). Observations were made by continuously 
scanning the telescope boresight in azimuth and elevation 
in a modified Lissajous pattern centered on 06h58m29.40s, 
-55d56m42.0s. The modified Lissajous pattern is defined as 
a function of time, t, by 

5Az = 5.5' sin at + 2' sin at/30. 
<5E1 = 5.5'sinfa + 2'sinto/30 

where a/6 = 8/9 and <5Az and oEl are physical coordinates 
relative to the field center. The actual values of a and 6 are 
normalized to limit the peak telescope velocity to 300"per 
second. A total of 61 maps of the cluster were completed, 
each taking 41.7 minutes, and collectively yielding an equiv- 
alent total on-source integration time of 193 s per 3" x 3" 
pixel. The data are reduced and co-added using the stan- 
dard AzTEC data analysis p ipeline with very similar tech- 
niques to those described in IScott et all l|2008h . The resul- 
tant map has an average point-source flux-error of 0.54 mjy 
over the inner 150 arcmin 2 area. A full analysis of the ex- 
tended 1.1 mm Sunyaev-Zel'dovich emission of the cluster 
and other point sources in the map will be presented in sub- 
sequent papers. 
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Figure 1. AzTEC image of MMJ065837-5557.0. The colorbar 
shows the shading in units of signal-to-noise (S/N) of a point 
source detection. The underlying noise in the map varies by 4.6% 
across the portion of the image shown. The faint diffuse emission 
to the west of the point source is dominated by the extended SZE 
signal from the Bullet cluster (Ezawa et al. in preparation) how- 
ever due to the point source filtering of the map, no morphological 
information about the SZE should be taken from this figure. 
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2.1 Pointing Corrections and Source Position 

A small correction to the telescope pointing model is ap- 
plied to all observations of the Bullet Cluster field based on 
periodic observations of a bright point source, PKS 0537- 
441, with a 1.1 mm flux of ~5 Jy. We derived the cor- 
rection by linearly interpolating the pointing offsets mea- 
sured from 4' x 4' maps of the point source observed every 
2 hours. Pointing measurements always bracketed observa- 
tions of the cluster. Our resulting pointing uncertainty, as 
measured from stacking the 1.1 mm flux at the locations 
of 212 known radio sources in our map of the GOOD S-S 
field (using the technique described in lScott et al.l l|2008h ) is 
4.8" rms. This is a random pointing error which broadens 
the PSF of the instrument by less than 0.5"and suppresses 
the point source response by an amount that is negligible 
compared to the photometric error. 

To identify the correct optical/lR counterpart to 
MMJ065837-5557.0 we next check for any systematic off- 
set in the pointing, which could be mechanical in nature or 
could be due to environmental influences on the telescope 
at the time of the observation. In either case, systematic off- 
sets are likely to vary from field to field. Without a known 
mm-bright object in our Bullet Cluster map to use as a 
pointing reference (or a sufficient number of fainter objects 
spread throughout the field to allow a stacking analysis) 
there is no way to determine any astrometric offset specific 
to this field. Instead, we estimate our overall astrometric un- 
certainty as the variance in residual pointing offsets (after 
the corrections applied from our observations of PKS 0537- 
441) measured from similar maps of 5 high redshift radio 
galaxies, along with the offsets measured from the stacking 
analysis from the GOODS-S field. Figure [5] shows the mea- 
sured offsets in right ascension and declination (physical co- 
ordinates) from these data. The thick red cross at ARA, 
ADec= (0.10", 0.18") shows the mean offset and la errors 
(3.24" in RA and 2.61" in Dec), i.e., consistent with no sys- 
tematic specific offset. Hence, in the analysis and discussion 
that follows, we apply no offset correction to the derived 
centroid of the point source. We assume that the standard 
deviation of systematic pointing offsets measured in other 
AzTEC maps, taken in a similar manner to the Bullet Clus- 
ter data, is representative of any positional offsets. 

To get the most accurate centroid position for 
MMJ065837-5557.0, we fit a 2-dimensional gaussian to the 
unfiltered image of the source. In addition to the source 
position we fit for an amplitude, the semi-major and semi- 
minor axes and position angle of the ellipse, and an arbitrary 
DC offset. Fitted parameters and their statistical errors are 
given in Table [T] 

Table 1. Derived parameters for MMJ065837-5557.0. The abso- 
lute calibration-error of the flux-density is given in parentheses. 



Flux density 

SZE-corrected flux density 
RA centroid position 
DEC centroid position 
Source FWHM (RA) 
Source FWHM (Dec) 
Position angle of elongation 
reduced x 2 1 °f fit 



15.9 ±0.5 (±1.3) mjy 
13.5 ±0.5 (±1.0) mjy 
06:58:37.31 ±0.02s 
-55:57:01.5 ±0.32" 
36 ± 1.3 arcseconds 
32 ± 1.2 arcseconds 
34 ± 8 degrees 
0.94 
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Figure 2. Measured positional offsets determined from maps of 
five high redshift radio galaxies (squares) and by stacking radio- 
detected AzTEC sources in the GOODS-S field (x). Error bars 
represent the la uncertainty in the centroid of the 1.1 mm flux of 
the radio galaxies. The red thick cross denotes the mean offset of 
the fields with bar lengths equal to the square root of the weighted 
variance of the set of radio galaxy offsets. 



2.2 Calibration and Flux Determination 

The AzTEC Bullet Cluster data are calibrated using Uranus 
as a prim ary calibrato r and adopting the technique de- 
scribed in IWilson et al.l (2008). Beammaps of Uranus were 
made twice each night. A linearly interpolated calibration 
factor is applied to Bullet Cluster observations taken be- 
tween Uranus beammaps. The calibration factor derived 
from the nearest Uranus beammap in time is used for Bullet 
Cluster observations taken before the first Uranus beammap 
of the evening and for Bullet Cluster observations taken 
following the last Uranus beammap of the evening. Since 
most of the science observations took place following the 
last Uranus beammap of the evening, we estimate an upper 
limit to our overall calibration uncertainty from the standard 
deviation in reported fluxes from 31 pointing observations 
of PKS 0537-441, which is 6% PKS 0537-411 is k nown to 
be a strongly varying source dRomero et all Il994h and so 
this upper limit should be considered robust. Adding this in 
quadrature with the 5% uncert ainty in the brightnes s tem- 
perature of Uranus at 1.1 mm (|Griffm fc Ortonlll993r ) gives 
a total calibration error of 8%. 

We use the 2-dimensional gaussian fit described above 
to estimate the flux density of MMJ065837-5557.0 in the 
absence of any foreground SZE emission from the cluster. 
Since the AzTEC map has been Wiener filtered to opti- 
mally identify point sources, an unknown degree of attenua- 
tion has been applied to the extended SZE signal. Thus the 
fitted point-source flux-density, allowing for a lower-limit on 
the contribution of SZE, is <Si.i mm = 13.5 ± 0.5(±1.0) mjy 
where the first error is the statistical error and the second 
(in parentheses) is the error in absolute calibration. 

The fitted size of MMJ065837-5557.0 is larger and 
slightly elongated compared to the beam size measured by 
performing the same fit to the 31 observations of the point- 
source PKS0537-441 (FWHM = 29 ± 1.2 arcseconds), sug- 
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gesting that the flux of MMJ065837-5557.0 may be due to a 
close pair, or possibly triplet (see §3.5), of confused sources. 
Until we have performed a more detailed subtraction of the 
SZE emission, or undertaken higher-resolution observations, 
we cannot rule out that the flux from MMJ065837-5557.0 is 
from a single point-source. 



3 DISCUSSION 

A comparison with the measured surface-densi ty of SMGs 
detect e d towards wi d e-area blank-fields (e.g. IScott et al.l 
(2008); ICoppin et alj l|2006h and references therein) is suf- 
ficient to demonstrate how unexpected is the detection of 
such a bright source. Although the submillimeter source- 
counts are unconstrained at these brightness-levels, a wide- 
area survey of S> 1 sq. degree is necessary before such an 
intrinsically-bright object could be expected to be discov- 
ered in a random survey. Two other similarly bright sources 
have been found in previous SCUBA surveys: J02399-0136, 
with an observed 850/im flux-density of 26 ± 3 mjy due 
to the am plification by the massive cluster Abell 370 at 
z — 0.37 (Ivi son et al. |l998); and the high-redshift radio 
galaxy 8C1909+722 which shows extended submillimetre 
emission and has an integrated flux of 34.9±3 mjy at 850/im 
l|Stevens et al.ll2003h . 

Given that MMJ065837-5557.0 has been detected near 
the peak in the weak-lensing map associated with the more 
massive of the two merging clusters (see Figure [3]), and that 
the Bullet Cluster shows optical signatures (i.e. arcs) of 
strong lensing, magnification of this source and amplifica- 
tion of the (sub)millimeter fluxes is a natural explanation of 
the high apparent brightness. Due to the low spatial resolu- 
tion of these single-dish ASTE observations (30" FWHM at 
1.1 mm) and the astrometric uncertainty, however, we must 
consider and dismiss other alternatives for the source of the 
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Figur e 3. The weak-lensing map (green contours), from 
I Clowe et alj ll2006f ) of the Bullet Cluster overlaid on an optical 
image that shows the two merging clusters of galaxies, with the 
less-massive "bullet" sub-cluster to the west of the main cluster. 
The white contours towards the peaks of the mass distribution 
indicate their positional uncertainty. The position of the bright 
AzTEC source MMJ065837-5557.0, which lies close to the peak 
of the more massive cluster, is shown as a white cross. The white 
circle indicates the scale of the AzTEC beam with 30" diameter 
FWHM. The full AzTEC map is significantly larger than the area 
of the optical image shown here. 



millimeter emission including Galactic stars, potential con- 
fusion with foreground Galactic cirrus, or galaxies intrinsic 
to the cluster at z — 0.297. 



3.1 Rejecting a Galactic origin for 
MMJ065837-5557.0 

Within a 6.5" radius 2a positional error-circle of the AzTEC 
source there is a compact point-source at the resolution of 
HST (object S in Figured]). The object is identified as a 
star in the Naval Observatory Merged Astrometric Dataset 
(NOMAD). S has an i?-band magnitude of 18.8, which corre- 
sponds to a flux-density of 0.19 mjy. Assuming a blackbody 
photospheric spectrum (T e s ~ 4500K) and that there is no 
extended circumstellar dust envelope, which is consistent 
with the declining flux-density of object S from 3.6 to 8 /im 
in the corresponding Spitzer imaging (see Fig. 3) , then S can- 
not be responsible for the 13.5 mjy emission at 1.1 mm. For 
completeness, we considered the possibility that the "stellar- 
object" S is an uncatalogued quasar at an unknown redshift. 
In the absence of an optical spectrum, we derive the IR col- 
ors of source S: S8.0/S4.5 < 0.40, S5.8/S3.6 ~ 1-7. Source 
S is bluer than a typical IR AGN and falls outside of the 
region of the colo r-color diagram shared by SMGs and AGN 
lYun et al] (|2008T l. We therefore reject source S as the source 
of the millimeter flux. 

At a Galactic latitude of b ~ —21° the foreground dust 
emission associated with Galactic cirrus is moderately-weak 
in the vicinity of the Bullet Cluster, with a 100/xm surface 
brightness of ~ 5 MJy/sr. Adopting the measured 60/100/im 
cirrus color-temperature of ~ 20K (and a dust emissivity- 
index f3 = 2), we estimate that the contribution to the 
measured AzTEC 1.1mm flux-density is ~ 1 mjy in the 
30" FWHM beam. This estimate is sensitive to the exact 
temperature of the cirrus, and can increase to 5 mjy for 
colder cirrus (T = 15K). However, the spatial- filter applied 
during the AzTEC data-reduction suppresses any extended 
foreground-emission on scales much larger than the point 
source response and therefore acts as a local sky-subtraction 
at the position of MMJ065837-5557.0. Unless the measured 
power due to cirrus is dominated by emission on angular- 
scales smaller than 30", foreground Galactic cirrus is un- 
likely to be responsible for the bright point-source emission 
of MMJ065837-5557.0. For the reasons given above we con- 
fidently reject a Galactic origin for MMJ065837-5557.0. 



3.2 Rejecting a cluster-member origin for 
MMJ065837-5557.0 

The optical galaxy nearest to the AzTEC centroid is an ellip- 
tical galaxy, which is a spectroscopically confirmed member 
of the Bullet Cluster. Located only 1.5" from the AzTEC 
source l ocation, this B = 21.58 elliptical galaxy at z = 
0.2958 |Barrena et all [20021 is a natural candidate for a 
counterpart to the 1.1mm continuum emission. It is also the 
nearest and the brightest IRAC source located within the 
la positional uncertainty of MMJ065837-5557.0. 

As shown in Figure [5] both the optical and the IRAC 
near-IR data of the galaxy are completely consistent with 
the spectral energy distribution (SED) of a typical elliptical 
galaxy, dominated by a predominantly old stellar popula- 
tion. The predicted 1.1mm flux density for such an elliptical 
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Figure 4. Left panel: HST ACS image (606W filter) showing a field of 23" X 23" centered on the centroid position (shown as a red cross) 
of MMJ065837-5557.0. The dimension of the cross indicates the ±1<t positional uncertainty of 3 arcsecs for MMJ065837-5557.0, which 
corresponds to a physical distance of ~ 22 kpc for a galaxy at 2 < z < 3. The dashed-lines, shown only on the HST and IRAC 8.0/xm 
images, bisecting the two red galaxy systems (A/C & B/D), and close to the position of the elliptical galaxy (E - a spectroscopically- 
confirmed member of the Bullet Cluster), are the approximate locations of the critical-line of infinite magnification for a galaxy at z = 2.7 
(adapted from Gonzalez et al. 2008). A stellar-object, S, is also identified. The optical galaxies A and B, enclosed within rectangles, 
are two candidates for highly-lensed background sources that may be high-,z counterparts to the AzTEC source MMJ065837-5557.0. 
The IRAC 3.6/xm (middle panel) and 8.0/xm (right panel) images of the same field as the HST image are shown with optical contours 
overlaid. Clear color differences and spatial displacement s , give n in Table 1, are seen between the optical components A and B, and the 
IR components C and D (shown as circles). iBradac et alj (|2006) present an alternative mass model in which the critical-line passes closer 
to object A. It requires higher-resolution (sub) millimeter observations with ALMA to confirm which of these components (if any) is the 
counterpart to MMJ065837-5557.0. Optical-IR spectroscopy are required to show whether the galaxy systems A/C, and B/D are multiple 
lensed-imagcs of the same background galaxy with heavy patchy dust-obscuration, or if components A,B,C and D are separate galaxies 
involved in strong-interactions that stimulate an episode of luminous starformation, or alternatively if they are physically-unconnected 
galaxies in our line-of-sight. 



galaxy, however, is nearly 4 orders of magnitudes smaller 
than that measured by AzTEC. However, major mergers of 
two gas-rich galaxies that lead to the ULIRG phenomenon 
are also known to pro duce an elliptical-like stellar remnant 
jHibbard fc Yunll 19991) . We explore the possibility that this 
elliptical galaxy may host a substantial and compact reser- 
voir of cold gas and dust by plotting the SED of Arp 220 (a 
prototypical ULIRG) scaled to match the measured 1.1mm 
flux density of MMJ065837-5557.0. The discrepancy be- 
tween the observed flux density and the scaled Arp 220 SED 
is only a factor of a few in the optical bands, but the dif- 
ference grows quickly and substantially at A > 1 fj,m. The 
absence of the PAH feature in the IRAC 8 /im band can be 
interpreted as evidence that the elliptical galaxy does not 
contain a large reservoir of cold gas and dust. This galaxy 
is not detected by the Infrared Astronomy Satellite (IRAS), 
and we show the 3<r upper limits of 120 mjy and 480 mjy in 
the IRAS 60 and 100 nm band|[] in Figure El An Arp 220- 
like SED that can account for the measured 1.1mm flux 
density can be ruled out with > 10cr significance. Therefore, 
we confidently reject the possibility that the elliptical galaxy 
is the origin of the extreme 1.1 mm continuum flux density 
of MMJ065837-5557.0. 

We can not confidently rule out other, more exotic al- 
ternative possibilities as the source of the bright millimeter 
emission such as nearby cooling-flow galaxies or exception- 
ally strong line emitters. As described below we instead focus 
on the most plausible explanation : two red IRAC sources 
(identified by iBradac et al.l ([2006) and further studied by 

1 obtained using ADDSCAN/SCANPI process avail- 
able through NASA/IPAC Infrared Science Archive 
(http://irsa.ipac.caltech.edu). 



lGonzalez~e t al. (2008])) representing a background galaxy or 
galaxies magnified by the Bullet cluster lens. 



3.3 MMJ065837-5557.0 - a highly-magnified, 
high-redshift starburst galaxy 

The Bullet Cluster is a high- velocity merger (~ 4500kms _1 ) 
of two massive systems (the main-cluster and the "bullet" 
sub-cluster) of galaxies taking place in the plane of the sky. 
Although there is continued discussion regarding the exact 
mass rati o (Mmain/Mb uiiet) of the two components ranging 
from 3:1 (jNusserlbOOsT ) to 10:1 (|Barrena et al.ll2002l ). the to- 
tal dark-matter hal o cluster mass is constrained to the rang e 
~ 0.8-2xl0 15 M Q (|Vikhlinin et alj|2006l ; lc"lowe et alj|2006l ). 
MMJ065837-5557.0 is within 15" of a peak in the underlying 
mass-distribution of the main galaxy-cluster and so a plau- 
sible scenario is that we are detecting amplified emission of 
a background sou rce close to this positi on. F urthermore, the 
lensin g models of lMehlert et alj l|200ll ) and lGonzalez et al.l 
(2008) show that critical- lines (caustics) of infinite magnifi- 
cation for galaxies at z — 3.24 and z = 2.7 respectively pass 
within a few arcsecs of the centroid-position of MMJ065837- 
5557.0 (Fig(4]). Thus, strong magnification of a background 
submillimeter galaxy is our preferred explanation for the 
unusual brightness of MMJ065837-5557.0. We now turn to 
the archival HST and IRAC data to identify the most likely 
counterparts of the mm-wavelength emission. 

Figure U shows a close-up view of the HST and IRAC 
fields near MMJ065837-5557.0, along with the approximate 
position of the lensing caustic-line of infinite magnification. 
Having ruled out the elliptical galaxy closest to the AzTEC 
centroid as the source of the strong millimeter wavelength 
emission, we focus on the two reddest IRAC sources (la- 
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Figure 5. Spectral Energy Distribution of the z = 0.297 cluster 
member galaxy E (FigHJ near the MMJ065837-5557.0 position 
centroid. The measured photometry data for the galaxy E are 
shown as filled squares, and t hey are compared wit h three ellipti- 
cal galaxy SED templates bv lPolletta et al] | |2007|) with different 
ages (2, 5, & 13 Gyr old shown as dotted, dashed and solid-lines 
respectively). The two 3cr upper limits from the IRAS 60 /an and 
100 /an bands are also shown. The Arp 220 SED by Polletta et al. 
is normalized to the 1.1mm AzTEC measurement (filled circle) 
for comparison. 



beled C and D in the IRAC 8 /im image of FigHJ) that lie 
within the uncertainty of the AzTEC centroid position, and 
that are also bisected by the caustic-line. The IRAC col- 
ors of objects C and D (see Table [2} are similar to those 
of dusty, optically-obscured, high-redshift SMGs. Objects C 
an d D also appear to be the same two galaxies identified 
by iBradac etatl i|200rj ) as candidates for galaxies at z > 6 
on the basis of their IR colors. We argue in § 13.41 that ob- 
jects C and D lie at significantly lo wer-redshifts, z ~ 2.7, a 
conclusion independently reached bv lGonzalez et al] (2008). 

A comparison of the HST image at 606 nm and Spitzer 
images at 3.6 and 8.0 /jm shows that two optical sources of 
different morphologies (objects A and B in Figure|4)| lie ^ 2" 
to the south of the IR centroids of objects C and D, which are 
most clearly defined in the 8.0pim image. A small positional 
offset of RA=— 0.4 arcsecs, Dec=0.0 arcsecs was applied to 
the archival HST data to align it with the IRAC image. 
The error in this offset, determined from the registration of 
stars detected in both fields, is insufficient to account for the 
apparent separation of the IR and optical pairs of galaxies. 
Before considering further the possibility of the association 
of the IR sources C and D with the optical HST sources A 
and B, and collectively all four sources with MMJ065837- 
5557.0, we present an estimate of the photometric-redshifts 
of the two IR sources. 



3.4 IR photometry, colors, and 
photometric-redshifts 

The difference in resolution between optical HST and IR 
Spitzer observations is significant. While the two extended 
red Spitzer sources (C & D) are clearly seen in the original 
8 fim IRAC images, the halo of the elliptical galaxy (object 
E in Figure [4} makes it non-trivial to measure their IR col- 
ors. Furthermore it is difficult to unambiguously determine 
if the IR emission of source C is associated with the optical 
emission of A, and similarly the IR source D with that of 
the optical source B. To minimize the contamination of the 
elliptical galaxy on the extracted fluxes of sources A, B, C 
and D, we model the elliptical galaxy using a 2-dimensional 
light-profile fit to the high-resolution HST image at 606 nm. 
A subtraction of this model from the corresponding 775 nm 
data did not reveal any significant residuals at the position 
of the elliptical. This same model was then smoothed and 
scaled to the resolution and peak intensity of the ellipti- 
cal galaxy emission in all 4 Spitzer IRAC bands and sub- 
tracted from each IRAC image, producing smooth residuals 
at the center of the elliptical that matched the surround- 
ing sky values. The optical fluxes at 606 nm and 775 nm 
for sources A and B were obtained using 1.5 and a 1.0 arc- 
sec radius aperture photometry, respectively, and adopting 
sky values estimated from concentric 0.5 arcsec-width rings, 
after masking any prominent objects, like star S. The IR 
fluxes for sources C and D at 3.6, 4.5, 5.8 and 8.0/im were 
extracted via aperture photometry after subtraction of the 
elliptical galaxy, as explained above. Star S was also mod- 
eled by a PSF, scaled and subtracted from the images to 
produce smooth profiles of the underlying object D, with- 
out creating large residuals in its structure. Sky values were 
estimated from areas without contamination of faint objects 
within concentric 3-8 pixel rings multiple times, and 2.5-3.0 
pixel (3.0-3.6 arcsec) radial aperture photometry was per- 
formed for each estimated sky value, and corrected to ac- 
count for flux losses using the PSF model corrections of the 
IRAC Data Handbook version 3.0. No attempt was made 
to formally measure the errors introduced in the photom- 
etry due to the elliptical and star subtractions, since the 
limited resolution of the IRAC images is not sufficient to 
constrain sub-pixel positional accuracies. Since both C and 
D lie outside the centroid of E and S, these additional errors 
are estimated to be less than 10%. All derived fluxes are 
listed in Tabled 

Figure[6]shows an IRAC color-color plot, S's.O/jm/^.s^m 
vs. S 5 .8um/S 3 .6um, fa IRAC-detected SMGs in the GOODS- 
N |Pope et al]|2006l ), Subaru-XMM Deep Field (Clements et 
al. 2008), and Lockman Hole (Dye et al. 2008) fields, along 
with the location of sources C and D. Both sources have IR 
colors consistent with other known SMGs. 

A crude photometric redshift for MMJ065837-5557.0 
can be estimated, based on the IRAC photometry, using 
an empirical polynomial-fit to the colors of similar galax- 
ies. This technique has been used successfull y in the optical 
and IR for various gal a xy po pulations (e.g. 1 Connolly et al] 
(|l995h and lPope etafl (j2006h ). We adopt a functional form 
z = a + E bi log Si , where a and hi are coefficients that are fit 
to the appropriate galaxy population (SMGs in this case), 
and Si are the IRAC flux-densities of the training galaxies at 
3.6, 4.5, 5.8 and 8.0^im. We fit our model to the IR colors of 
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the 16 SMGs with unambiguous single optical/IR counter- 
parts, complete IRAC photometry, and robust spectroscopic 
redshifts in the GOODS-N and SHADES fields (Pope et al. 
2006; Dye et al. 2008; Clements et al. 2008). A best fit results 
in the relationship: 

z = 1.6393 + 3.9134 log(S 5 .8/S 3 .6)- 2.3490 log(5 , 8 .o/5 , 4.5)(l) 

In order to test the accuracy of this photometric- 
redshift estimate for the SMG population we use a leave- 
one-out cross-validation technique, where each member of 
the training-galaxy set is systematically excluded from the 
fit, and then that best-fit solution is used to derive its photo- 
metric redshift. This ensures that the accuracy of this pho- 
tometric redshift technique is determined independently of 
the galaxies used to derive it. Figure [7] shows the compari- 
son of the spectroscopic redshifts of the training set to the 
photometric redshifts derived from the IRAC colors. We es- 
timate the error in each photometric redshift by bootstrap- 
ping on the photometric errors in the colors. We find an av- 
erage error between photometric and spectroscopic redshifts 
of < 5z/(l + z) >= 0.15 for this technique. Although the 
accuracy of this method decreases as the redshift increases, 
at least based on this small sample, there is no apparent bias 
to it. 

Using this technique we find that both IR sources, C 
and D, have photometric redshifts of 2.7 ±0.2 where the un- 
certainty is derived from the observed colours of the galax- 
ies and bootstrapping on their individual errorbars. This 
supports the arguments that they are either multiply-lensed 
images of the same source, or that they are two galaxies 
separated (in the image-plane) by ~ 35 kpc and are under- 
going an interaction that stimulates star-formation in dust- 
obscured regions. 



3.5 Counterparts to MMJ065837-5557.0 

We now explore a number of scenarios regarding the possible 
associations of the optical and IR sources A, B, C and D 
(Fig.©. 

Scenario 1: Gradients of patchy dust obscura- 
tion (A=C and B=D). If we assume that the optical HST 
and IR Spitzer emission of components A and C, and sim- 
ilarly for components B and D, arises from different and 
distinct regions of the same galaxy, and that the galaxies 
A/C and B/D are both at redshifts z ~ 2.7 (SCO), then 
the physical separations of the IR and optically-emitting 
regions in each galaxy are < 16 kpc in the plane of the 
lens. We consider how this apparent separation decreases in 
the source plane. Using previously pu blished mass-models 
derived from strong lensing data alone (|Mehlert et al.|[200lh 
along with combined strong and weak-lensing re construction 
for the cluster 1E0657-558 (jBradac et al J 120061 ). the image 
plane is mapped back t o the source plan e using the inversion 
software LENSTOOL (| Julio et alj|2007h . We model the bul- 
let cluster as two large-scale smooth potentials, and include 
the elliptical galaxy, a confirmed cluster-mem ber, as an ad- 
dition al small-scale local perturber. Following iBradac et al.l 
(2006), the main cluster is modelled with a surface mass- 
density (or covergence) n(r) oc r -1,2 and velocity dispersion 
o~ ~ 1400 kms -1 , whilst the sub-cluster has a slightly shal- 
lower slope of k[v) oc r -0,9 and a ~ 1200 kms -1 . Typically 
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Figure 6. Color-color p lot of IRAC-dete cted submm galax- 
ies in the GOODS -N jPope et al] l2006h . SHADES SXDF 
Elements et alj|200Sl) and Lockman-Hole fields tove et al]|200Sl) 
with secure spectroscopic redshifts, represented as diamonds, 
where the hue of the symbols is proportional to the redshifts. 
SXDF850.21, at SkcWS^m = 8.27, Ss.s^/Sa.e^m = 1.33, 
•Zspec = 0.044 is out of bounds of the represented plot, but has 
been included in the analysis. The squares represent the two IR 
bright components (C and D) that are possible counterparts at 
z ~ 2.7 to MMJ065837-5557.0. The lines represent the color 
changes for three model spectral energy distributions of dusty 
galaxies (Siebenmorgen & Kriigel 2007) placed at increasingly 
larger redshifts, as indicated by the small labels. 
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Figure 7. Photometric vs. spectroscopic redshift of submm 
galaxies in the GOODS-N and SHADES fields. The photometric 
redshift for each galaxy has been calculated using an empirical 
fit to the colors of SMGs, following the method described in §3.3. 
The dispersion of values between spectroscopic and photometric 
redshifts is < 6z/(l + z) >= 0.15. 

the local amplification, k 2 — 7 2 , where 7 is the shear, pro- 
vides an comparable expansion of the image-plane compared 
to the source-plane. In the vicinity of the MMJ065837-5557.0 
however the shear-field of the cluster is perpendicular to the 
elongation of the cluster. Hence, adopting the above com- 
posite mass-model, the ~ 2arcsec separation of the pairs of 
optical and IR galaxies (A/C and B/D) in the image-plane 
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Table 2. Optical and IR photometry of galaxies within the positional uncertainty of MMJ0659-557. Sources A-E refer to those 
objects identified in Figure [4] All fluxes are given in ^tJy. Redshifts with superscripts p denote photometric-redshift estimates (this 
paper), and with s denotes a spectroscopic-redshift from Barrena et al. (2002). 



source 


redshift 


606nm 


775nm 


3.6/xm 


4.5/^m 


5.8^tm 


8.0/xm 


notes 


A 




0.8 ±0.05 


1.0 ±0.03 












B 




0.6 ±0.05 


1.3 ±0.03 












C 


2.7 ±0.2? 






27 ± 1 


44 ±1 


57 ±7 


54 ±2 


after E subtraction 


D 


2.7 ±0.2? 






25 ± 2 


41 ± 1 


80 ±8 


100 ±2 


after E subtraction 


E 


0.297 s 


63 ±4 


80 ±4 


192 ± 12 


135 ± 9 


112 ±5 


75 ±4 





is reduced only by a factor of 0.5 - 0.75, or a separation of 
1.0 — 1.5 arcsec (10 — 13 kpc), in the source-plane. 

It is therefore possible that the galaxies identified in the 
HST image are the same as those detected by Spitzer IRAC, 
i.e. A=C and B=D, implying that the thermal emission at 
millimeter wavelengths arises from distant, dusty starburst 
galaxies with a gradient of heavy, but patchy obscuration. In 
support of this argum ent recent K-band ima ging of the same 
composite systems bv lGonzalez et all (|2008i ) shows that the 
centroid of the 2.2/im emission lies between the optical HST 
and longer-wavelength 3-8/rni Spitzer data. 

Scenario 2: Interacting pairs of red &: blue galax- 
ies (A^C and B^D.) Alternatively the HST and Spitzer 
data may be indicating that we have two close pairs of 
strongly-interacting or merging-galaxies, in which each pair 
contains an optical galaxy and an IR dust-obscured star- 
burst. The observed separation of sources A and C (and 
equivalently B and D), i.e. ~ 8 arcsecs, corresponds to 
^ 30 kpc in source plane. The association of composite 
(possibly interacting) systems, red galaxies and blue com- 
panions, with the environments of SMGs is not uncommon 
|lvison et al . 1 120021 ) . The Antennae galaxy provides an ex- 
cellent local example of the complex wavelength-dependent 
morphology in a merging-galaxy with spatially-distinct re- 
gions of luminous optical emission from young stars, and 
even greater luminosity at millimeter wavelengths origi- 
nating from a young starburst population that remains 
heavil y-embedded in the du s ty, gas-ric h ISM (see for ex- 
ample Iwhitmore et all (| 19991 ) ; IXu et all (|2000l ); IWang et ail 
(12004 )1. 

Scenario 3: Two multiple lensed-images of a 
background galaxy (A=B=C=D). The final scenario 
that we consider, which is consistent with the cluster lens 
model, is one in which the pairs of optical and IR galaxies 
are bisected by a critical-line of infinite magnification. In 
this case objects A, B, C, and D may all be multiple lensed- 
images of the s ame background s ource . Similar examples 
are presented bv lKneib et al.l (|2004t ) and lBorvs et al l (|2004l ) 
in the clusters Abell 2218 (z = 0.17) and MS0451.6-0305 
(z = 0.55) respectively. In these cases bright submillime- 
ter sources detected by SCUBA (9-17 mjy at 850/im) are 
associated with multiple images of highly-magnified, multi- 
component optical and IR galaxies which are bisected by 
the critical- lines. After demagnification, the intrinsic fluxes 
of these particular SMGs are 0.4-0.8 mjy at 850^m, making 
them some of the faintest SMGs detected to date. 

An argument against the optical sources A and B be- 
ing multiple images of the same background galaxy is that 
they have different morphologies in the high-resolution HST 



observations (Figure [!}. Source A shows more extended and 
elongated emission than source B, which appears to con- 
sist of two compact emitting-areas. However, the presence 
of the cluster-member elliptical galaxy (source E in Fig[4j| 
is able to provide sufficient gravitational mass to distort the 
resolved optical images and create additional magnification 
of the background source (e.g. Natarajan & Kneib 1997). 
It is not necessary therefore that the two lensed images of 
the background source should have identical magnifications 
and morphologies. Given the optical brightness of galaxies A 
and B, this scenario can be easily tested via a spectroscopic- 
determination of their redshifts. 

Finally, we mention that our gravitional-lens model also 
predicts a third counterpart to the multiply-imaged back- 
ground source with only a modest magnification factor of 
3-5 at 06h 58m 33.5s, -55d 57m 29 s (J2000). A faint red 
IRAC source is seen at this position with no optical counter- 
part. We do not discuss the possibility of this third source 
further, except to mention that it is marginally consistent 
with the orientation of the slightly-extended profile of the 
millimeter emission from MMJ065837-5557.0. 

Regardless of the degree of association between the 
IR and optical sources, the IRAC photometric-redshifts are 
well- within the redshift distribution of the blank-field SMC 
population, and thus we believe that the most plausible ex- 
planation for the mm-wavelength flux is that 1) the AzTEC 
source MMJ065837-5557.0 lies behind the Bullet Cluster, 
2) it is lensed by the Bullet Cluster, and 3) it is a prod- 
uct of the combined emission of the IRAC objects C and 
D. This scenario also naturally explains the approximate 
equal distance from the AzTEC centroid position to the two 
IR galaxies. In the absence of a detailed lensing model in 
the vicinity of the AzTEC source and/or spectroscopic con- 
firmation, we assume in the remainder of this paper that 
both pairs of optical/IR sources (A/C and B/D) are associ- 
ated with MMJ065837-5557.0, which is a lensed background 
SMG. 



3.6 Contribution of the Sunyaev-Zel'dovich effect 
to the observed flux of MMJ065837-5557.0 

Two classes of effects can serve to increase the measured 
brightness of detected point sources: multiplicative effects 
(lensing) and additive contributions to the flux due to mm- 
wavelength emission along the line of sight to the source. 
MMJ065837-5557.0 clearly suffers from both. In this section 
we describe the non-negligible millimeter-wavelength emis- 
sion due to the thermal Sunyaev -Zel'dovich Effect (SZE, 
ISunvaev fc Zeldovichl (|l970l . Il972l )) that is added to the 
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intrinsic lensed emission of MMJ065837-5557.0. The ther- 
mal SZE is produced by inverse-Compton scattering of the 
CMB photons off electrons in the hot and dense inter-cluster 
medium (ICM) of clusters. We do not consider the contri- 
bution from the weaker kinetic SZE that is due to the bulk 
motion (velocity) of the cluster plasma with respect to the 
CMB reference- frame. Consequently the thermal SZE is a 
dominant source of secondary-fluctuations on the surface 
brightness of the CMB on small angular scales, less than 
a few arcminutes. 

Unable to perform a multi-wavelength spectral decom- 
position of the SZE and the bright background galaxy 
MMJ065837-5557.0, with AzTEC observations only at 
1.1mm, we can estimate the level of the thermal SZE con- 
tamination in two independent ways. 

First, taking advantage of the spatial-resolution of the 
AzTEC observations and assuming that the SZE emission 
from the cluster is smooth on spatial scales smaller than that 
of the AzTEC beam (30"), we can spatially-filter the AzTEC 
map and fit for the 1.1mm flux consistent with only point- 
source emission (as discussed in §2.1). Alternatively, as de- 
scribed below, we can model the SZE emission from the Bul- 
let cluster with a simple model for the electron-temperature 
and density distribution to estimate the 1.1mm flux from 
the SZE at the position of the AzTEC point-source. 

The derivation of the intensity of the SZE and its 
observed properties at radio to millimeter- wavelengths 
have be en reviewe d exten sive ly llSunvaev fc Zeldovichl 
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1999 
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2007). Briefly, considering only the 



thermal distortion, the frequency-dependent SZE is given 
by the simple expression 



SIv = Io\g{x)y), 



(2) 



where I = 2(k B T ) 3 /(he) 2 and T = 2.725K is the mean 
temperature of the CMB. The Comptonization parameter, 
y, equal to the optical-depth times the fractional energy-gain 
of each scattering event, is 
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n e T e dl 



(3) 



Finally the dimensionless spectral function, g(x), where 
x — hv/kBTo, describes the frequency-dependence of the 
thermal SZE. Neglecting relativistic corrections, g(x) has 
the form 



9{x) = 



e x + 1 



I) 2 



(4) 



Since the intensity of SZE is dependent only on the line- 
of-sight integral of the electron pressure through the Bullet 
cluster, we describe below our estimations for the local value 
of the electron temperature and density at the position of 
MMJ065837-5557.0. 

It is common to adopt an isothermal assumption within 
the cluster density distribution (typically a f3- model), or 
at least a weak radial-dependence on the tem perature of 
the in tracluster-medium within the virial radius (|lrwin et alj 
1999). Support for th is assumption is also found in hydro dy- 
namical simulations l|Eke et a l. 1998; Pearce ct al. 20q3). In 
the case of the Bullet cluster, Mar kevitch et al. I (|2002l ) have 
provided a temperature-map that shows a gas-temperature 



variation across the merging clusters between 10-24 keV, ig- 
noring the concentrated cold-spot or "bullet" (with a tem- 
perature of 6-7 keV). The spatial resolution of AzTEC (30" 
FWHM beamsize) is sufficient to estimate that the local 
temperature of the X-ray-emitting gas, at the position of 
MMJ065837-5557 0, is 14-16 keV (regions f-g in Fig.2 of 
IMarkevitch et all l|2002l )'). or similar to the average cluster 
temperature of ~ 15 keV. 



Whilst the density-profile of the gas to the west of the 
main X-ray peak (in the direction of the "bullet") requires 
the inclusion of a shock-front and an additional compo- 
nent to represent the lower-mass merging cluster, the sit- 
uation is less-complicated to the east of the X-ray peak 
brightness, in the direction of MMJ065837-5557.0. The lo- 
cal electron-density, n e , in the eastern hemisphere (at the 
position of MMJ065837-5557.0) is calculated using a single- 
component spherically-symmetric /3-model (/3=0.7) centered 
at the primary peak of the observed X-ray surface-brightness 
distribution. We begin by considering a rang e for the to- 
tal c l uster-mass, M to t = -8 - 1.9 x 10 15 M m ijClowe et all 
120061 ; IVikhlinin et al.ll2006T ). Following iKav et all |200ll ) we 
then assume that the ICM is fully ionized, with a helium 
mass fraction of 0.24, such that n e — 0.88p/m,H , where 
m,H is the mass of a hydrogen atom. The baryonic mass 
density, p, is calculated ass uming a global baryon fraction 
f b = M b /Mtot = 0.06/i" 3/2 (|Ettori fe Fabianll 19991 ). The ex- 
ercise was then repeated assuming an alternative /3-model 
centered on the peak of the weak-lensing signal from the 
main cluster that dominates the total cluster mass (§.3.3). 
Finally, provided that the electron temperature, T e , is simi- 
lar to the X-ray gas temperature, we determine that the con- 
tribution of the SZE at 1.1mm at the position of the AzTEC 
point-source MMJ065837-5557.0 is ~ 1.5 - 4.2 mjy/beam. 
This estimated range of underlying flux due to the SZE is 
consistent with what we measure directly from the extended 
component in the spatially-resolved 1.1mm emission at the 
position of the AzTEC point-source (§2.1). 



We therefore conclude that the 1.1mm flux-density of 
MMJ065837-5557.0 due solely to continuum emission from 
a lensed high-redshift SMC is in the range 12-14 mjy. An 
analysis of the detected extended AzTEC emission towards 
the Bullet Cluster (Ezawa et al. in preparation) will provide 
more accurate information on the location of MMJ065837- 
5557.0 with respect to the peak of the SZ effect, and conse- 
quently a more accurate measurement of the point-source 
contamination to the integrated SZE (a potential prob- 
lem for lower-resolution experiments) than the estimate de- 
scribed above. 



Assuming a range of spectral-indicies a ~ 3 — 4 (where 
F„ <x i/ a ) and dust temperatures to describe the rest- 
frame spectral energy distribution of MMJ065837-5557.0 in 
the millimeter regime, we also determine that MMJ065837- 
5557.0 is too faint to be detected in the low-resolution 
(FWHM~ 4.5') ACBAR obser vations at 1.4 mm given the 
published noise-level of 30/iK (I Gomez et alj|2004h . This is 
consistent with the lack of comment by Gomez et al. about 
a detectable signal at the null-frequency of the SZE. 
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3.7 The intrinsic luminosity of 

MMJ065837-5557.0 - a dusty high-redshift 
LIRG 

The magnification map of the Bullet Cluster by 
iMehlert et al.l (|200lT > suggests a robust and conservative 
lower-limit of A = 20 for the amplification of MMJ065837- 
5557.0. More typical amplifications of A = 1 — 3 have 
been found for SMGs behind low-z massive clusters 
jLedlow et all |2002| ; iKnudsen et al.l 120061 ). It is the prox- 
imity of MMJ065837-5557.0 to a lensing caustic-line that 
makes this such an exceptional case. The implication for 
such a large magnification is that the AzTEC source has 
an intrinsic flux of Si,i mm < 0.68 mjy, after subtraction of 
the expected contribution from the Sunyaev-Zel'dovich ef- 
fect. This is a factor of > 3 lower than the flux of the faintest 
blank-field SMGs detected in the full 225 sq. arcmin AzTEC 
field towards the Bullet Cluster (which occupies only the 
central 12 sq. arcmins). 

Assuming that the AzTEC point-source is a background 
galaxy at z ~ 2.7, then then intrinsic millimeter-wavelength 
brightness of MMJ065837-5557.0, after the corrections for 
the SZE contribution and amplification, corresponds to a 
FIR luminosity of < 1O 12 L0. Given the relative insensitiv- 
ity of the millimeter-wavelength flux-density as a function 
of redshift, this is a robust result provided that the AzTEC 
source lies in the redshift range ~ 1 < z < 8. Hence, 
rather than having discovered an extremely rare exam- 
ple of a hyperluminous dust-obscured galaxy, MMJ065837- 
5557.0 is more typical of a luminous IR galaxy (LIRG) with 
Z/fir ^ 10 11 — 1O 12 L0 in the high-redshift Universe, with 
a star-formation rate (SFR) in the range 5 — 50M Q /yr _1 
ijKennicutdl 19981 ). 

The conservative lower-limit on the magnification im- 
plies that even a more modest SFR may be appropri- 
ate, closer to that of only a mildly-active starburst galaxy. 
This would support the scenario in which the millime- 
ter emission arises from a distant starburst with patchy 
dust obscuration and a low dust-content, which in turn 
increases the likelihood of detecting both the optical and 
IR emission from the source (i.e. A/C and/or B/D). Only 
high-resolution (sub)millimeter interferometric imaging of 
MMJ065837-5557.0, to unambiguously identify the counter- 
part (s), combined with spectroscopic measurements to de- 
termine the redshifts (and possibly velocity dispersion) of all 
optical and IR components (A, B, C, D) will allow progress 
in understanding the nature of this potentially rare exam- 
ple of a low to intermediate luminosity millimeter-selected 
galaxy in the high-redshift Universe. 



4 CONCLUSIONS 

Sensitive continuum observations at 1.1mm with the AzTEC 
camera installed on the 10-m diameter Atacama Submillime- 
ter Telescope Experiment have been made towards the Bul- 
let cluster. An extremely bright point-source (MMJ065837- 
5557.0) is detected with an observed flux-density at 1.1mm 
of 15.9 mjy. The centroid position of MMJ065837-5557.0 lies 
close t o the largest mass- density peak in the weak-lensing 
map of lClowe et all <|2006h . 

Archival IR (Spitzer IRAC) images show the presence 



of two red sources, separated by ~ 8 arcsecs (or ^ 30kpc 
in the source-plane), that lie within the positional-error 
of MMJ065837-5557.0, and that also straddle the critical 
m agnification-lin e (i.e. a region of extreme amplifications) 
of IMehlert et al.l (|200ll ). The Spitzer IRAC colors of these 
two sources are typical of high-z SMGs, and photometric- 
redshifts place both sources at z ~ 2.7 ± 0.2. The proximity 
of HST-identified optical sources to both of these Spitzer 
IRAC galaxies leaves open 2 possibilities: 1) that each opti- 
cal/IR pair are individual galaxies undergoing close interac- 
tion, or 2) that each pair is in fact only one galaxy with a 
complex morphology of patchy dust-obscuration. 

Having rejected a Galactic origin and a cluster mem- 
ber origin as the source of the AzTEC emission, the above 
conditions suggest that MMJ065837-5557.0 is a background 
source associated with one or both of the red IR Spitzer 
sources, and is strongly-lensed by this massive cluster. The 
agreement in th e location of the crit ic al-line in the magni - 
ficat ion maps of IMehlert et al.l l|200ll ). iBradac et ail (|2006T) 
and i Gonzalez et al l Jgxri T^ Cluster place a strong 
conservative limit of > 20 for the magnification of the mil- 
limeter emission from MMJ065837-5557.0. There exist many 
similarities between this example of an extremely bright mil- 
limeter source, and the possibility of gravitationally-lensed 
interacting galaxie s that produce a b right submillimeter arc 
in MS0451.6-0305 jBorvs et al.ll2004h . 

After subtracting the contribution (^ 25%) due to the 
extended SZE at 1.1mm from the total AzTEC continuum 
flux of MMJ065837-5557.0, we estimate that the intrinsic 
(i.e. demagnified) FIR luminosity of MMJ065837-5557.0 is 
< 10 12 Lq provided the source lies in the redshift range 
~ 1 < 2 < 8. Hence MMJ065837-5557.0 is representative of 
the fainter population of starforming galaxies that is unde- 
tectable, without the assistance of lensing, with the current 
generation of small single-dish (sub)millimeter telescopes 
due to source-confusion and lack of sensitivity. Until the 
commissioning of the next-generation of larger millimeter- 
wavelength facilities (e.g. LMT, ALMA) the only opportu- 
nity to study this fainter luminosity class of IR galaxies, that 
dominate the production of the cosmic infrared background, 
is to take advantage of the significant amplification that ex- 
ist towards the critical magnification-lines of massive nearby 
clusters. 
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